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Center  in  the  understanding  of  water-entry  phenomena.  The  research  reported  herein 
was  supported  entirely  by  NAVSEA  Code  63R31_ under  cask  SR0230100.  The  author 
would  like  to  acknowledge  Dr.  Thomas  Pierce  of  NAVSEA  for  his  advice  and  interest 
in  this  program. 
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INTRODUCTION 


The  development  of  naval  weapons  requiring  delivery  through  an  air-water 
interface  necessitates  an  understanding  of  the  various  water-entry  phenomena  which 
have  a  predominant  influence  on  the  trajectory  of  a  water-entry  vehicle  from  initial 
impact  until  the  attrition  of  the  water-entry  cavity.  To  address  this  need  the 
large  Hydroballistics  Facility  and  Associated  Pilot  Tank  were  constructed  at 
the  Naval  Surface  Weapons  Center/White  Oak  Laboratory  specifically  for  the  study 
of  water-entry  phenomena.  A  program  was  initiated  in  1967  to  obtain  basic  vertical 
water-impact  acceleration  data  for  conical  shapes.  This  effort  was  later  extended 
to  include  ogives,  cusps  and  the  sphere.  A  digital  program,  Reference  1,  has 
recently  been  written  to  predict  the  forces  during  water-entry.  This  code  requires 
the  use  of  approximate  procedures  for  estimating  the  wetting  factor  (height  of 
the  heaved  water  surface  wetting  the  cone)  and  for  estimating  the  shape  of  the 
cavity. 

A  concurrent  program  was  also  initiated  in  1967  to  study  the  behavior  of 
the  water-entry  cavity.  In  a  recent  study,  Reference  2,  an  analytical  model  was 
developed  to  predict  the  dynamic  behavior  and  shape  of  the  water-entry  cavity, 
neglecting  the  effect  of  the  free  surface.  This  model  successfully  predicts  the 
frequency  response  of  the  cavity  which  occurs  after  cavity  closure.  The  model 
was  developed  from  energy  considerations  and  assumed  that  the  fluid  velocity  was 
in  a  direction  normal  to  the  velocity  of  the  water-entry  vehicle.  Near  the  free 
surface  this  is  obviously  not  true.  The  compliance  of  the  free  surface  results 
in  a  large  vertical  velocity  component  as  evinced  by  the  heaving  of  the  free 
surface  and  the  splash  that  occurs  in  the  neighborhood  surrounding  the  point  of 
water  impact.  In  extending  the  study  to  include  the  effect  of  the  free  surface, 
a  theoretical  model  for  predicting  the  heaving  motion  of  the  free  surface  was 
developed.  This  model  successfully  predicts  the  wetting  factors  for  cones  obtained 
experimentally  as  reported  in  Reference  3.  The  use  of  this  model  in  the  digital 
program  of  Reference  1  should  eliminate  the  need  for  the  approximate  methods 
presently  employed. 


^Wardlaw,  A.  B. ,  Morrison,  A.  M.  and  Baldwin,  J.  L. ,  "Prediction  Of  Impact 
Pressures,  Forces,  and  Moments  During  Vertical  And  Oblique  Water-Entry",  NSWC/ 
WOL/TR  77-16,  January  1977. 

2 

Smith,  C.  W. ,  "Experimental  Investigation  of  the  Behavior  of  Vertical  Water-Entry 
Cavities",  NSWC/WOL/TR  78-135,  September  1978. 

^Baldwin,  J.  L. ,  "Vertical  Water-Entry  of  Cones"  NOLTR  71-25,  February  1971. 
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The  ability  of  the  model  to  predict  the  wetting  factors  for  cones  encouraged 
an  extension  of  the  model  to  include  a  method  for  predicting  the  impact  drag  force 
during  vertical  water-entry.  Successful  correlation  of  the  cone  data  reported 
in  Reference  3  was  obtained  after  incorporating  a  momentum  correction  to  the 
theoretical  model.  Only  two  empirical  constants  were  necessary  to  correlate  the  cone 
data  reported  in  Reference  3  (10,  15,  20,  30,  45,  60,  90,  120  and  140  degree  cones). 
The  model  was  then  extended  to  incorporate  a  changing  prismatic  coefficient  (for 
displaced  volume).  Experimental  data  for  the  sphere,  Reference  4,  and  for  ogives 
and  cusps,  Reference  5,  were  then  correlated.  For  the  sphere  and  each  ogive  and 
cusp,  a  different  set  of  empirical  constants  was  required.  The  author's  hope 
that  one  set  of  empirical  constants  might  correlate  all  geometries  was  not  realized. 
Correlation  of  the  empirical  constants  appears  to  be  a  possibility,  however,  test 
data  for  additional  geometries  must  be  obtained  and  correlated  before  this  can 
reasonably  be  attempted. 


See  footnote  3  on  page  5 
a 

Baldwin,  J.  L.  and  Steves,  H.  K. ,  "Vertical  Water-Entry  of  Spheres",  NSWC/WOL/ 
TR  75-49,  May  1975. 

^Baldwin,  J.  L. ,  "Vertical  Water-Entry  of  Some  Ogives,  Cones,  and  Cusps",  NSWC/ 
WOL/TR  75-20,  March  1975. 
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THEORETICAL  MODEL  FOR  THE  HEAVED  FREE  SURFACE 


Shiffman  and  Spencer,  in  Reference  6,  employed  an  approximating  ellipsoid 
to  obtain  an  expression  for  the  wetting  factor  and  drag  coefficient  for  the  vertical 
water-entry  of  cones.  Their  model  assumed  that  the  approximating  ellipsoid  should 
have  a  fineness  ratio,  Yi  such  chac: 

Y  =  i  =  cot  (j)  (1) 

This  assumption  implies  that  the  velocity  of  the  approximating  ellipsoid  is  equal 
to  che  velocity  of  penecration  of  the  cone.  In  modifying  Shiffman  and  Spencer's 
analysis  to  obtain  better  agreement  with  the  experimental  data  for  cones  in 
Reference  3,  no  assumption  was  made  for  the  fineness  ratio  of  the  ellipsoid. 

Figure  1  defines  che  nomenclature  used  to  describe  the  geometry  of  the  cone  and 
ellipsoid.  In  this  figure,  che  depth  of  penecration  of  the  ellipsoid  is  defined: 

A  =  a  -  (h-H)  (2) 

and  the  velocity  of  penecration  is  obtained  by  differentiating  equation  (2): 

A  =  a  -  (h-H)  (3) 

To  account  for  Che  compliance  of  che  free  surface  ic  was  assumed,  as  did 
Shiffman  and  Spencer,  chac  che  heaved  free  surface: 

n  (r,  c)  »  -  z  (4) 

is  locaced  in  che  equatorial  plane  of  che  ellipsoid  (x  =  o).  In  this  plane 
che  velocity  of  che  fluid,  and  therefore  che  free  surface,  is  vertically 
upward.  On  che  intersection  of  che  heaved  free  surface  wich  che  ellipsoid, 

Che  fluid  velocity,  from  equation  (A-29)  of  Appendix  A,  is: 

h  -  H  =  f|  |  =  K  A  (5) 

r=b 

where  is  defined  in  Appendix  A,  equation  (A-20) . 

A  third  equation  is  obtained  by  performing  a  mass  balance; 


JO 

2irp  J"  qrdr  +  piTb^  (h-H)  =  pirb^  C  h 


( o ) 


where  it  is  assumed  that  the  correlation  factor,  CD,  is  the  prismatic  coefficient 
of  the  entering  body. 
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Shiffman,  M.  and  Spencer,  D.  C.,  "Force  of  Impace  On  a  Cone  Striking  a  Water 
Surface  (Vertical  Entry)",  Comm.  Pure  and  Appl.  Mach.,  Vol  IV  (1951). 
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Equation  (6)  can  be  written,  after  substituting  equation  (A-35)  from  Appendix  A  for 
the  first  term: 


A  +  (h-H)  «  C  fi 
P 


(7) 


Since,  for  an  arbitrary  configuration,  b  and  Cp  can  be  written  as  functions 
of  h,  it  is  more  convenient  to  use  h  as  the  independent  variable.  Equations  (3), 
(5)  and  (7)  then  become,  respectively: 

a'  =  a'  -  (1  -  h'  )  (8) 


1  -  H  -  K  A'  (9) 

A  +  ( 1  -  H  )  -  C  (10) 

P 


To  obtain  an  expression  for  a  and  y,  A  from  equation  (8)  is  substituted  into 
equation  (10): 

a'  .  C  (11) 

P 

integrating: 

h 

a  =•  f  C  dh  (12) 

P 

o  r 

and  dividing  by  b: 

y  ■t-f  /  CP  dh  !13> 

o 

To  obtain  an  expression  for  h-H  and  the  wetting  factor,  equacion  (9)  is 
integrated: 


h  -  H  =  /  K.  A  dh 


Substituting  equation  (9)  for  (1  -  H  )  into  equation  (10): 


(14) 


(1  +  K.)  A  *  C 
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wnere  b.0  is  defined  in  Appendix  A,  equation  (A-22),  as  a  function  of  the  fineness 
ratio,  y.  Equation  (iw)  can  now  be  written: 


h  -  H  = 


16) 


and  the  wetting  factor: 


w  = 


h 

H 


1 


C  dh 
P 


For  a  cone,  C 

P 


1/3,  b  =  h  tan  (|) 


and  equation  (13)  becomes: 


1  h  _  1  _  ,Q, 

Y  “  3  b  _  3  COC  (2 


(17) 


(18) 


For  a  given  cone  angle  a,  y  and  c_  are  constants  and  equation  (17)  becomes: 
1 

w  =  t - rr 

1  -  Q  /6 
o 


(19) 


The  calculated  wetting  factors,  using  equation  (19),  are  tabulated  in  Table  1  and 
compared  with  the  experimental  data  from  Reference  3.  The  predictions  are  in 
good  agreement  with  the  experimental  results  and  well  within  the  variation  of 
experimental  results. 

MODEL  FOR  THE  IMPACT  DRAG  FORCE 


To  obtain  the  impact  drag  force,  the  principle  of  conservation  of  momentum 
is  employed.  Due  to  rotational  symmetry,  the  total  momentum  of  the  fluid  is 
directed  vertically  downward.  From  equation  (A-40)  of  Appendix  A: 


P 


MH  *  2ltp / JlJrdrdx 


a  Kx  A 


(20) 


An  expression  for  the  'added  mass'  is  obtained  from  equation  (20): 


M 


aX 


where : 


K 


1  A 

i  ■  ■ 

H 


(21) 


(22) 


The  force  on  the  body  entering  the  water  is  directed  vertically  upward  and 
is  equal  to: 

I 

•  •  *  M  *2 

F-P-MH-4L.H  (23) 

H 

3 

See  footnote  3  on  page  3 
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TABLE  1 

WETTING  FACTORS  FOR  CONES,  w  *  h/H 


Cone 

Experimental  Results, 

Ref.  3 

Predicted 

Values 

Angle,  a 

Average 

Variation,  X 

Eq . ( 19) 

(w-w)/w 

Degrees 

w 

Lowest 

(-) 

Highes  t 

(+) 

w 

% 

10 

- 

1.0275 

- 

15 

- 

1.0465 

- 

20 

1 . 0609 

10.39 

10.34 

1.0655 

+  .43 

30 

1.1038 

10.33 

11.42 

1.1018 

-  .18 

45 

1.1290 

11.25 

7.71 

1.1507 

+1.92 

60 

1.1661 

8.65 

11.73 

1.1937 

+2.37 

90 

1 . 2406 

8.19 

13. 11 

1.2687 

+2.27 

120 

1.3101 

11.45 

11.42 

1.3381 

+2.14 

140 

1.3977 

4.34 

13.35 

1.3858 

-  .85 

f 

t 


] 


1 


w  ***** nwwtft, 


'll  'i  no*—— P»i—H  i 
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where  n  is  assumed  chat  the  velocity,  H,  is  a  constant.  The  instantaneous  impact 
drag  coefficient  is  defined: 


Cd  = 


2F 


2M 


4  (b“aXj 

2  •?  2  '  ~  3  2  r 

pitr  H  p  TTr  H  r  H 

o  o  o 


(24) 


For  a  cone,  X  and  y  are  constants,  b  =  h  tan  (j) ,  and  equation  (24)  can  be 
written : 


C  .ii,  <iL)2 
d  3  '  r 

H  o 


The  maximum  impact  force  occurs  when  b  =  r 


(25) 


max 


4  X 
1 


(26) 


H 


The  maximum  impact  drag  coefficients  reported  in  Reference  3  for  cones  could 
not  be  adequately  correlated  with  the  use  of  equation  (26).  To  obtain  a  model  for 
correlating  impact  drag  data,  a  momentum  term  was  added  to  equation  (20).  In 
Figure  1  the  expanding  geometry  of  the  water-entry  body  was  neglected.  To 
accommodate  the  increasing  radius  of  the  water-entry  body,  a  radial  velocity 
field  proportional  to  b  is  required.  To  create  this  local  radial  velocity  field 
it  was  assumed  that  an  augmented  vertical  velocity  field  proportional  to: 


U  *  *\/(b)2  +  (h-H)2 


(27) 


is  deflected  by  an  additional  impact  force.  As  shown  in  Figure  2,  the  change 
in  the  vertical  velocity  field  is  proportional  to: 


AU  *  *\/(b) 2  +  (h-H)2  -  (h-H) 

and  the  total  momentum  of  the  fluid  is  now  assumed  to  be: 
P  *=  Ax  MH  +  A2  MAU  “  MH  [ai  +  A2  (^H)  ] 


(28) 


(29) 


where  Aj  and  A2  are  empirical  constants  to  be  obtained  through  the  correlation 
of  experimental  data.  The  total  impact  force  on  the  water-entry  body  becomes: 


F  -  P 


P  -  MH  [  Aj  +  A2  r]  +  A2  MH  f 


h2  {m'  [  +  a2  r]  +  a2  m  r'}  /h' 


where : 


r  =  c-^)  -  {V<b’>2  +  a  -  h')2  -  a  -  h’)}  /h 

H 


(30) 


(31) 


3 

See  footnote  3  on  page  5. 
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FIGURE  2  VELOCITY  FIELD  DEFLECTION 


.3 


NSWC  TR  80-243 


14 


NSWC  TR  80-243 


where : 


(39) 

(19) 

(40) 

(18) 


The  cese  data  in  Reference  3  for  the  impact  drag  force  exhibited  this  parabolic 
behavior.  It  was  necessary,  therefore,  to  correlate  only  the  maximum  drag 
coefficients: 


C 


d 


max 


(41) 


The  results  of  the  correlation,  with  A.j  =  .447  and  A2  *  2.1052,  are  tabulated 
and  compared  with  experimental  data  in  Table  2.  The  correlation  is  in  good 
agreement  with  the  experimental  results  and  is  well  within  the  variation  of  the 
test  data. 


Test  daca  for  the  sphere,  obtained  from  Reference  4,  were  correlated.  Two 
correlations  were  obtained  and  tabulated  in  Table  3  for  comparison  with  experimental 
results.  Both  correlations,  the  first  with  A3  =  0,  and  che  second  with  A3  =  A2 . 
are  in  good  agreement  with  the  test  data  (generally  within  2.2%).  The  daca  were 
correlated,  using  least-square-fit,  for  penetration  depths  up  to  H  =  .40  r0. 

At  this  penetracion  depth  the  heaved  surface  on  the  sphere  is  located  at  h  =  .658 
rQ  where  che  slope  on  Che  sphere  is  20  degrees.  Experimental  daca,  according  to 
Hsu-Perry  and  reported  in  Reference  7,  indicates  that  flow  separation  occurs  in 
this  neighborhood.  When  separation  occurs,  che  impact  drag  force  coefficients 
are  no  longer  applicable.  The  drag  force  goes  through  a  transition  phase,  from 
impact  drag  to  cavicy- running  drag. 

The  method  used  to  define  che  geometry  of  ogives  and  cusps  is  shown  in  Figure 
3,  and  is  given  che  name:  a/8  ogive.  The  hemisphere  can  be  defined  as  a  90/90 
ogive  and  the  cone  can  be  defined  as  a  a/o  ogive. 


See  footnote  3  on  page  5 
See  footnote  4  on  page  6 

I 

May,  A.,  "Wacer-Entry  and  che  Cavity-Running  Behavior  of  Missiles1',  SEA  HAG/ 
TR  75-2,  1975. 
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TABLE  2 

MAXIMUM  IMPACT  DRAG  COEFFICIENTS  FOR  CONES,  Cd 

max 

Ai  *  .4470 

Cone  Experimental  Results,  Ref.  3  Correlation,  A2  m  2.1052 


Angle,  a 

Average 

Avariation,  % 

Degrees 

5d 

max 

Lowest 

(-) 

Highest 

(+) 

10 

.0221 

19.5 

19.1 

15 

.0427 

9.3 

9.3 

20 

.0696 

11.1 

14.5 

30 

.1277 

5.2 

7.7 

45 

.2858 

8.6 

7.1 

60 

.4608 

8.5 

9.4 

90 

1.175 

5.3 

4.9 

120 

2.840 

7.2 

9.7 

140 

5.547 

4.8 

3.9 

Eq.  (41)  (Cd  -  Cd  )/Cd 

max  max  max 


cd  % 


max 

.0221 

+  .4 

.0425 

-  .4 

.0677 

-2.7 

.1331 

4.2 

.2733 

-4.4 

.4766 

+3.5 

1.187 

+1.0 

2.815 

-  .9 

5.483 

-1.2 

IMPACT  DRAG  COEFFICIENT  FOR  SPHERE,  Cd 


FIGURE  3  GEOMETRY  OF  OGIVES  AND  CUSPS 
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The  cest  daca  for  ogives  and  cusps,  obcained  from  Reference  5,  were  correlated 
and  tabulated  in  Tables  4  thru  7  for  the  88.4/32,  60/43,  88.84/-31  and  59.94/-43.4 
ogives.  The  correlations  with  A3  =  0  and  A3  ■  A2  appear  to  be  comparable.  For 
the  ogives,  both  correlations  are  in  good  agreement  with  the  test  daca  and, 
excluding  one  daca  point,  are  well  within  the  variation  of  the  test  daca.  The 
correlations  for  the  cusps,  although  not  always  within  the  variation  of  the 
experimental  daca,  appear  to  be  reasonably  good. 

The  empirical  constants  obcained  from  the  leasc-square-fit  of  che  impact 
drag  force  data  are  tabulated  in  Table  8.  The  empirical  constants  for  A3  =  0 
are  also  plotted  in  Figure  4.  Correlation  of  Che  empirical  constants  appears 
to  be  a  possibility,  however,  test  data  for  additional  geometries  must  be  obtained 
and  correlated  before  this  can  reasonably  be  attempted. 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  successful  correlation  of  Che  experimental  data  for  che  impact  drag 
coefficients  for  cones,  ogives,  cusps  and  che  sphere  indicates  the  usefulness 
of  Che  drag  model  for  correlating  the  impact  drag  coefficients  for  a  wide  range 
of  water-encry  configurations.  The  data  base  for  water-entry  forces  should  be 
extended  and  correlated.  Hopefully,  this  would  lead  to  che  successful  correlation 
of  Che  empirical  constants.  The  usefullness  of  having  a  general  correlation  model 
for  predicting  che  impact  drag  forces  during  water-entry  is  obvious.  Impact  drag 
forces,  for  which  chere  is  no  cest  daca,  could  then  be  predicted  with  reasonable 
confidence  in  the  results.  A  method  for  predicting  the  drag  force  in  the  transition 
phase,  from  impact  drag  to  cavity- running  drag,  should  also  be  developed.  A  drag 
model  for  che  cavity-running  phase  has  previously  been  developed  and  reported 
in  Reference  8. 


5 


See  footnote  5  on  page 


6 


g 

Gaula,  J.  R.  and  Smith,  C.  W.,  "A  Drag  Model  For  Cavitacing  Axisymmetric  Bodies", 
NOLTR  73-173,  Sept  1,  1973. 
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IMPACT  DRAG  COEFFICIENT  FOR  60/43  OGIVE 
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TABLE  8 

SUMMARY  OF  EMPIRICAL  CONSTANTS 


Ogive 

A1 

A2 

a/6 

A3  =  0 

A3  =  A2/2 

a3  =  a2 

a3  «  0 

A3  =  A2/2 

A3  *  A, 

59.94/-43.4 

2.2334 

2.5897 

2.8006 

2.6238 

1 . 9346 

1.5305 

60/0 

.447 

.447 

.447 

2.1052 

2.1052 

2.1052 

60/43 

.3967 

.7309 

.9984 

2.1656 

1.9252 

1.7328 

88.84/-31 

2.1726 

2.7841 

3.1134 

2.2136 

1.5389 

1.1781 

90/0 

.447 

.447 

.447 

2.1052 

2.1052 

2.1052 

88.4/32 

.2729 

.6785 

1.0282 

2.2174 

2.0536 

1.9123 

90/90 

.6254 

.9111 

1.2619 

.8820 

.9724 

1.0828 
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APPENDIX  A 

POTENTIAL  FLOW  SOLUTION  FOR  THE  AXIAL  MOTION  OF 
ELLIPSOIDS  OF  REVOLUTION 


Semi-elliptic  coordinates  (p,  £ ,  oj)  are  defined  in  Reference  9  in  terms  of 
the  cylindrical  coordinates  (x,  r,  u)): 

x  =  E  p  C  (A-l) 

r  =  E  C*/l  ~  U2  (A- 2) 

a;  =  U!  (A-3) 

Where  E  and  defined  in  Table  A-l,  are  introduced  to  alleviate  the  need  to 
distinquish  between  ovary  and  planetary  ellipsoids. 

The  surface  of  the  moving  ellipsoid,  X,  *  C0>  is  defined: 


<!>■ 


1 


( A- 1 0 ) 


Substituting  for  x  and  r  from  equations  (A-l)  and  (A-2) ,  the  following  values  for 
£  and  £  on  the  surface  of  the  ellipsoid  are  obtained: 


o 


(A-ll) 


*o 


(A-12) 


The  line  elements  in  semi-elliptic  coordinates  have  the  value,  Reference  9: 


(A-13) 


9 

Munk,  M.  M.,  "Ellipsoids  Of  Revolution",  Chapter  VII,  Division  C,  Volume  I,  Of 
"Aerodynamic  Theory",  W.  F.  Durand,  Dover  Publications  Inc.,  N.  Y.,  1963. 
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TABLE  A-l 
NOMENCLATURE  FOR 
OVARY  AND  PLANETARY  ELLIPSOIDS 


Planetary,  b>a 

Ovary,  a>b 

E2 

+  (b2  -  a2) 

-  (b2  -  a2) 

(A-4) 

;2  ♦  1 

C2-l 

(A-5) 

V2 

2  2 

5  +  U 

c2  -  p2 

(A-6) 

f(?) 

cot  1  Z 

5  l"  C^T 

(A- 7) 

F(C) 

-[  f  -  2  «'>] 

(A-8) 

G(^) 

cos"1  Y 

ln[  y  +Vy2  ~  l] 

(A-9) 

V1  -  Y2 

■\/y2  -  1 
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ds 


C 


ds 


Ev  d£ 

l 


E  C  V  1  -  d^ 


( A- 1 4 ) 

(A-15) 


Where  V  is  defined  in  Table  A-l.  The  velocity  components,  in  terms  of  the  velocity 
potential,  $,  are: 


V 

y 


=  V 1  -  y2  3£ 
3S  Ev  3u 

y 

3£  _  4  3$ 

3S^  Ev  3; 


( A- 1 6 ) 


(A-l 7) 


r  -  -  -  _  1  3$ 

to ' "  %  ’ "  e  ** 


( A- 1 8 ) 


The  velocity  potential  for  the  axial  motion  of  an  ellipsoid  in  an  infinite 
inviscid  fluid  is  well  known,  Reference  9.  The  velocity  potential  can  be 
expressed: 


♦  *  E  Kj  A  p  ?  F(?)/F(?o) 


(A-19) 


Where  is  the  inertial  factor  for  axial  motion: 
K,  -  — 2- 


1  2  -  a 


(A-20) 


and: 


“o  "  50  V  FV 


(A-21) 


or,  in  terms  of  the  fineness  ratio,  y  m  ^  : 

“o  ■  rH?  [ 1 ' Y  e(*>  ] 


(A- 22) 


F(£)  and  G(,y)  are  defined  in  Table  A-l. 

Substituting  equation  (A-19),  for  the  velocity  potential,  into  equations  (A-16) 
thru  (A-18),  the  velocity  field  is  obtained: 


Ki  A 


f(; 


;  f(;) 


(A-23) 


See  footnote  9  on  page  26 
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V 

C 


Kx  A 


(A- 24) 


V 


U) 


0 


(A-25) 


1c  is  assumed  that  che  plane  of  the  heaved  free  surface: 

n  (r,  c)  *  -  z  ( A-26) 

is  in  che  equatorial  plane  of  che  moving  ellipsoid: 

U  *  o  (A-27) 

In  this  plane  che  velocicy  of  che  fluid,  and  therefore  che  free  surface,  is 
vertically  upward  (see  Figure  1): 

A  *  ■  -  Vy.  -  K,  A  ?(;)/F(Co)  (A-28) 

U=o 


On  che  intersection  of  che  heaved  free  surface  and  Che  ellipsoid,  Q  *  £  : 


Kx  A 


(A-29) 


The  mass  flow  race  ac  che  heaved  free  surface  can  be  expressed: 

00 

m  -  2irp  /  f|  rdr  (A-30) 

b 

An  expression  for  rdr  can  be  obtained  from  equation  (A-2) 
rdr  ■  E2  £ ( l  -  y2)  Cd£  -  £2  ydy  ^ 

On  che  surface,  y  3  o,  equation  ( A—  31)  becomes: 
rdr  *  E2  £d£ 

Substituting  equations  ( A— 3 2 )  and  (A-28)  into  equation  (A-30): 

•  00 

m  -  2up  E2  |k~  y  f  F(0  ^dQ  (A-33) 

“•o  c 
*o 


( A-  3  1 ) 

(A-32) 
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Since : 


;  r^j  ;d;  =  - 

5o 

ar.d  from  equacions  (A-21),  (A-20)  and  (A-12)  respectively  : 


(A-34) 


F(5  ) 

o 


^0 


2 


equation  (A-33)  becomes: 

GO 

ih  =  2up  /  fi  rdr  =  pirb^  A 
b 


(A-35) 


Due  to  rotational  symmetry,  the  total  momentum  of  the  fluid  is  directed 
vertically  downward: 


P  *  MH  =  2iip  JJ  V  rdr  dx  =  -  2vp  ff  |£  rdr  dx 

y=o  5-*oo 

*  2irp  /  d  ]  rdr  +  2np  J  <J>]  rdr 

v*=i  p“° 

-  2~p  /  rdr  (A-36) 

jc*«° 

On  the  surface  x  a  and  on  the  heaved  free  surface,  y  *  0,  the  velocity  potential 
is  zero  (<J>  «  0)  and  equation  (A-36)  becomes: 

.  y=o 

P  *  MH  ■  2irp  /  4>  I  rdr  (A-37) 

y-1  «e 

On  the  surface  of  the  ellipsoid,  ?  *  £o>  equations  (A-19)  and  (A-31)  become, 
respectively: 


♦  I  '  E  (  K  i  y  '  a  K,  A  W 


(A-38) 
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